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ABSTRACT 

Uncertainties regarding the phylogenetic history within the Bivalvia have impeded attempts to 
understand evolution within the group. Estimating the evolutionary relationships surrounding the 
Trigonioida has been especially problematic and has led to disparate hypotheses regarding (1) 
the origin and subsequent diversification of unionoid bivalves and (2) autobranch gill character 
state transitions. In order to test alternative hypotheses of trigonioid phylogeny, 613 base pairs 
of DNA sequence of the cytochrome c oxidase subunit I gene were analyzed from 14 species of 
bivalves (ingroup) and three species of non-bivalve mollusks (outgroup). All phylogenetic analy¬ 
ses, using either nucleotide or inferred amino acid sequences, produced trees that robustly 
placed Neotrigonia (Trigonioida) as the sister taxon to a monophyletic Unionoida. Furthermore, 
the Autobranchia, Veneroida, and Mytiloida were supported as monophyletic groups in these 
analyses, whereas the Bivalvia was not. These phylogenetic relationships suggest that (1) there 
was a single invasion of freshwater by a unionoid bivalve ancestor, (2) trigonioid rather than ven- 
eroid bivalves gave rise to the Unionoida, (3) either the eulamellibranchous or filibranchous gill 
condition has evolved multiple times within the Autobranchia, and (4) the molluscan bivalved 
body plan may have evolved more frequently than traditional phylogenetic hypotheses suggest. 

Key words: Bivalvia, Trigonioida, Unionoida. cytochrome c oxidase I, mtDNA, phylogenetics, 
convergence. 


INTRODUCTION 

The higher-level phylogenetic relationships 
within the Bivalvia are poorly understood at 
present (e.g., Allen, 1985). This statement is 
corroborated by the many disparate hypothe¬ 
ses of evolutionary relationships that have 
been proposed for the higher taxa within the 
Bivalvia (e.g.. Purchon, 1990; Waller, 1990; 
Starobogatov, 1992; Cope, 1996; B. Morton, 
1996; Salvini-Plawen & Steiner, 1996). This 
plethora of phylogenetic hypotheses likely 
stems, in part, from arbitrary choices to ex¬ 
clude certain types of characters and a gen¬ 
eral lack of explicit and rigorous phylogenetic 
analyses. Although major anatomical charac¬ 
ter suites distinguish two of the nominal sub¬ 
classes within the Bivalvia, that is, Proto- 
branchia and Autobranchia, the ordinal-level 
relationships within each subclass are not 
strongly supported by multiple shared-derived 
morphological features (e.g., Waller, 1990). 

Two contrasting phylogenies of the higher- 
level relationships within the Bivalvia have 


recently been proposed by Salvini-Plawen 
& Steiner (1996; Fig. 1A) and Waller (1990; 
Fig. IB). The former study is an exemplar in 
that explicit data matrices and phylogenetic 
methodologies were presented. A major point 
of disagreement between the two phylogenies 
presented in Figure 1 is the placement of the 
Trigonioida, a relatively ancient bivalve taxon 
that, although currently depauperate, was tax- 
onomically diverse during the Mesozoic (Cox 
et al., 1969: Allen, 1985). The Salvini-Plawen 
& Steiner hypothesis (Fig. 1A) indicates that 
the trigonioids are most closely related to pte- 
riomorph bivalves (represented herein by 
mytiloids), with the Veneroida closely related 
to the Unionoida. Their proposed sister taxon 
relationship for the trigonioid and pteriomorph 
bivalves was supported by the shared pres¬ 
ence of (1) byssate larvae and adults and (2) 
abdominal sense organs in these taxa 
(Salvini-Plawen & Steiner, 1996). Alterna¬ 
tively, the hypothesis of Waller (Fig. IB) indi¬ 
cates that the Trigonioida is the sister taxon to 
the Unionoida. However, no evidential sup- 
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FIG. 1. Bivalve relationships, based on morphological characteristics, according to Salvini-Plawen & Steiner 
(1996; A) and Waller (1990; B). Both hypotheses indicate a single origin for the Bivalvia (1) and filibranchous 
lamellibranch gills (2) but differ in their implications for eulamellibranch gill evolution. Hypothesis A indicates 
a single origin for eulamellibranch gills (3) while B indicates either two origins (solid bars) or a single origin 
followed by a reversal to the filibranch condition in trigonioids (hatched bars). 


port for this relationship was provided by 
Waller (1990). 

The evolutionary affinities of the Trigonioida 
have been a contentious subject for over a 
hundred years (e.g., Steinmann, 1888: Neu- 
mayr, 1889; Cox, 1960; B. Morton, 1987; 
Healy, 1989, 1996) and are central to an un¬ 
derstanding of bivalve mollusk character evo¬ 
lution, since representatives of the single ex¬ 
tant genus ( Neotrigonia) have a mixture of 
seemingly primitive (e.g., filibranchous gills 
[possessing ciliary-linked gill filaments], lack 
of posterior mantle fusion, nacreous shells; 
Cox, 1960; Allen, 1985: B. Morton, 1987) and 
derived features (e.g., multi-vesicular sperm 
acrosome; Healy, 1989). This mosaic of 
seemingly primitive and derived character 
states has contributed to the erection of multi¬ 


ple, disparate hypotheses of trigonioid evolu¬ 
tionary affinities. According to the hypothesis 
of Salvini-Plawen & Steiner (1996: Fig. 1A), 
the placement of the trigonioids is consistent 
with a single origin of eulamellibranch gills 
(tissue-linked gill filaments) from filibranch 
gills. In contrast, the phylogenetic placement 
of trigonioids in the hypothesis of Waller 
(1990; Fig. IB) suggests that there were ei¬ 
ther (1) two origins of eulamellibranch gills or 
(2) a single origin of eulamellibranch gills fol¬ 
lowed by a reversal to filibranch gills in the 
trigonioids. 

The phylogenetic uncertainty within the Bi¬ 
valvia impedes the rigorous testing of evolu¬ 
tionary hypotheses. Therefore, the establish¬ 
ment of a robust phylogenetic hypothesis for 
higher taxa within the Bivalvia will enable criti- 
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cal evaluations of (1) trigonioid evolutionary 
relationships and (2) hypotheses of auto¬ 
branch bivalve gill character state transitions. 
Molecular systematic analyses have been 
useful in situations where morphological 
analyses were inconclusive (e.g., Avise, 
1994). Therefore, mitochondrial DNA(mtDNA) 
sequences, obtained from the cytochrome c 
oxidase subunit I (COI) gene, were used to 
construct a phylogenetic hypothesis for the 
major bivalve lineages represented in Figure 
1. COI was chosen for this analysis because of 
its slow rate of evolution relative to other mito¬ 
chondrial protein coding genes, relative ease 
of unambiguous sequence alignment (e.g., 
Brown, 1985; Simon et al., 1994; Russo et al., 
1996), and demonstrated appropriateness for 
this particular analysis. 


MATERIALS AND METHODS 
Organisms 

The molluscan taxa used in this study, with 
GenBank accession numbers and primary lit¬ 
erature citation (where applicable), are as fol¬ 
lows: (1) ingroup, Class Bivalvia, Subclass 
Protobranchia, Order Solemyoida, Solemya 
velum (U56852), Order Nuculoida, Nucula 
tenuis (U56851), Subclass Autobranchia, 
Order Unionoida, Mutela rostrata (U56849), 
Amblema plicata (U56841), Anodonta cygnea 
(U 56842), Margaritifera margaritifera 
(U56847), Order Trigonioida, Neotrigonia 
margaritacea (U56850), Order Mytiloida, 
Modiolus modiolus (U56848), Geukensia de- 
missa (U56844), Order Veneroida (the five 
veneroid COI sequences are from Baldwin et 
al. [1996]), Corbicula fluminea (U47647), 
Rangia cuneata (U47652), Mercenaria mer- 
cenaria (U47648), Mytilopsis leucophaeata 
(U47649). Dreissena polymorpha (U47653), 
(2) outgroup. Class Scaphopoda Dentalium 
sp. (U56843), Class Gastropoda Lepetodrilus 
elevatus (U56846), Class Polyplacophora 
Katharina sp. (U56845). 

Methods 

Total DNA was isolated from somatic (man¬ 
tle) tissues of nine species of bivalves. Male 
gonadal tissues were specifically avoided dur¬ 
ing dissections to prevent comparisons of 
non-orthologous sequences (due to the actual 
or potential presence of doubly uniparental in¬ 
heritance of mtDNA in some bivalve taxa; for 


example. Skibinski et al., 1994; Zouros et al., 
1994; Hoeh et al., 1996). DNA was also iso¬ 
lated from representatives of three additional 
molluscan classes (i.e., Gastropoda, Polypla¬ 
cophora, and Scaphopoda) for use in gener¬ 
ating outgroup sequences. Subsequently, a 
710bp fragment of COI was PCR amplified 
and cycle sequenced for each of the 12 taxa 
as described elsewhere (Folmer et al., 1994). 
Both strands of the COI fragment were se¬ 
quenced from each of two individuals from 
each terminal taxon to guard against PCR- 
based contamination artifacts. The resulting 
12 COI sequences, plus the five veneroid bi¬ 
valve COI sequences from Baldwin et al. 
(1996), were readily aligned by eye using 
MacClade 3.05 (Maddison & Maddison, 
1992). Sixteen of the 17 OTUs produced se¬ 
quences 613 bp in length, while that of 
Geukensia was 616 bp. The increased length 
of the Geukensia sequence was due to an au- 
tapomorphic, three nucleotide (single codon) 
insertion event. These three contiguous nu¬ 
cleotides, which are phylogenetically uninfor¬ 
mative for the taxa considered herein, were 
deleted prior to all phylogenetic analyses. No 
additional hypothesized insertion or deletion 
events were necessary to obtain the align¬ 
ment utilized in the subsequent analyses. 

The suitability of the COI data set for phylo¬ 
genetic analyses at the required hierarchical 
level was evaluated by plotting the substitu¬ 
tion pattern of transitions and transversions 
for each codon position (e.g., Orti & Meyer, 
1996). Furthermore, the degree of phyloge¬ 
netic signal within the COI data set was eval¬ 
uated using the g 1 statistic of a random tree 
distribution (from 100,000 random trees; e.g., 
Hillis, 1991; Hillis & Huelsenbeck, 1992) as 
implemented in PAUP (Swofford, 1993). Phy¬ 
logenetic analyses were carried out on the 
COI nucleotide sequences using the maxi¬ 
mum likelihood ([ML]. DNAML in PHYLIP 
3.5c; Felsenstein, 1993), neighbor-joining 
([NJ], MEGA 1.02; Kumar et al., 1993), and 
maximum parsimony ([MP], PAUP 3.1.1; 
Swofford, 1993) algorithms. Katharina . the 
only non-conchiferan mollusk taxon, was 
used to root the resulting topologies. A transi¬ 
tion/transversion ratio of 2.0 was utilized in 
the ML analyses and the gamma distance 
(alpha = 0.5, using the Tamura-Nei model of 
nucleotide sequence evolution) was used to 
generate the pair-wise genetic distances for 
the NJ analyses. This particular distance 
takes into consideration among-site substitu¬ 
tion rate variation (e.g., Yang, 1996). Further- 
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more, MP and NJ (again using gamma dis¬ 
tances, alpha = 0.5) analyses were conducted 
on the inferred COI amino acid sequences 
(using the Drosophila mtDNA genetic code). 
Multiple random terminal taxon addition order 
runs, combined with global branch rearrange¬ 
ment options, were employed to generate 
topologies from ML and MP analyses. These 
options increased the probability of finding the 
actual best topology under each of the two op¬ 
timality criteria (e.g., Hendy et al., 1988; Mad- 
dison, 1991). The robustness of the resulting 
topologies was evaluated by bootstrap analy¬ 
ses (1,000 replicates for MP and NJ. 100 
replicates for ML). 

The best COI-based topology derived from 
the DNAML analysis was compared with the 
phylogenetic hypotheses presented in Figure 
1 using the Kishino-Hasegawa test (paired z 
test; Kishino & Hasegawa, 1989) as imple¬ 
mented in DNAML. To this end, the topological 
constraints option in PAUP was used to gen¬ 
erate 82 user trees (all trees < five steps longer 
than the shortest trees found by PAUP) repre¬ 
senting the two tree topologies (for the partic¬ 
ular taxa evaluated herein) in Figure 1. Each of 
these user trees was then compared to the 
best DNAMLtree by the Kishino-Hasegawa al¬ 
gorithm. This test evaluated the significance of 
any potential incongruence between the mor¬ 
phology- (Figure 1) and COI-based topologies 
(Swofford et al., 1996). In addition, character 
optimization, using MacClade (Maddison & 
Maddison, 1992), was carried out on the COI- 
based topologies to investigate their implica¬ 
tions for morphological character evolution 
within the Bivalvia. 


RESULTS 

Scatter plots of the relationship between 
the number of transitional and transversional 
substitutions at each of the three codon posi¬ 
tions and the proportion of nucleotide differ¬ 
ences (all positions) for the COI sequences 
revealed that only transitional substitutions at 
the third codon position had reached satura¬ 
tion (Fig. 2). Because saturated categories of 
substitution can contribute to erroneous esti¬ 
mates of evolutionary history (e.g., Swofford 
et al., 1996), all first and second position sub¬ 
stitutions together with only transversions at 
the third codon position were included in the 
COI nucleotide data matrix used for phyloge¬ 
netic analyses. Of the 613 nucleotide posi¬ 
tions in the transformed COI data matrix, 229 


were invariant, while 315 (102 from 1st codon 
positions; 59 from 2nd; 154 from 3rd) were 
phylogenetically informative using the parsi¬ 
mony criterion. Analysis of the tree length dis¬ 
tribution of 100,000 randomly generated 
trees, using all 17 taxa, suggests that there is 
a significant amount of hierarchical structure 
within the transformed COI data set (g 1 = 
-0.794; with 384 variable sites, p « 0.01; 
Hillis & Huelsenbeck, 1992). This analysis 
was repeated using only the Corbicula, Den- 
talium, Katharina , Lepetodrilus, Modiolus, 
Neotrigonia, and Solemya COI sequences 
(seven taxa) in order to minimize the number 
of closely related taxa present in the hierar¬ 
chical structure analysis. Significant structure 
was still present in this truncated COI data 
matrix (g 1 = - 0.527; with 295 variable sites, p 
< 0.01), which suggests that the hierarchical 
structure present in the original data matrix 
was not solely due to the presence of closely 
related taxa. The findings from the plots of 
substitution pattern and g 1 statistics are con¬ 
sistent with the hypothesis that significant 
phylogenetic signal exists in the transformed 
COI nucleotide data matrix and validate the 
latter’s use in this particular phylogenetic con¬ 
text (e.g., Swofford et al., 1996). Since amino 
acid substitution rates are lower than the un¬ 
derlying nucleotide substitution rates (e.g., Li 
& Graur. 1991), it follows that the inferred COI 
amino acid sequences are not saturated and, 
therefore, also appropriate for this level of 
phylogenetic analysis. 

The best tree produced by ML analysis of 
the transformed COI nucleotide matrix is pre¬ 
sented in Figure 3 along with bootstrap per¬ 
centages for the NJ (above branches, 1,000 
replicates), MP (below branches. 1,000 repli¬ 
cates), and ML (in parentheses, 100 repli¬ 
cates) analyses. Only bootstrap percentages 
greater than 50% are shown. This topology is 
largely congruent with the best trees pro¬ 
duced by MP (three equally parsimonious 
trees) and NJ analyses (trees not shown). 
However, the topological relationships of the 
non-autobranch taxa on the best MP trees 
were identical to those portrayed in the trees 
derived from analyses of the inferred COI 
amino acid sequences (see below). 

The Kishino-Hasegawa test results (Table 
1) indicate that the morphology-based topol¬ 
ogy represented in Figure lAwas significantly 
worse (p < 0.05) than the best topology from 
ML analyses of the transformed COI data ma¬ 
trix (Figure 3). However, the morphology- 
based topology represented in Figure 1B was 
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Proportion of Nucleotide Differences (all positions) 

FIG. 2. Substitution pattern for transitions (s) and transversions (v) at each codon position for the COI frag¬ 
ment analyzed herein. The number of transitions and transversions is plotted against the proportion of nu¬ 
cleotide differences over all positions. 


not significantly worse (p > 0.05) than that de¬ 
rived from ML analysis. These test results 
confirm significant discordance between the 
morphology-based topology represented in 
Figure 1A and the COI-based topology (e.g., 
Swofford et al.. 1996). 

Figure 4 represents the best topology gen¬ 
erated from NJ analysis of the inferred COI 
amino acid sequences with bootstrap per¬ 
centages for NJ (above branches) and MP 
(below branches) analyses. Only bootstrap 
percentages greater than 50% are shown. 
The strict consensus tree (not shown) of the 


39 equally parsimonious trees produced by 
MP analysis of the inferred COI amino acid 
matrix is less resolved but entirely congruent 
with the NJ tree. 

All analyses of the transformed COI nu¬ 
cleotide and inferred amino acid matrices 
were consistent with the monophyly of the Au- 
tobranchia, Veneroida, Mytiloida, and Union- 
oida (e.g., Figs. 3, 4). The trees in Figures 3 
and 4 are especially noteworthy due to the 
strong inference that the filibranchous 
Neotrigonia (Trigonioida) is more closely re¬ 
lated to eulamellibranchous freshwater mus- 
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FIG. 3. Best tree topology produced by maximum likelihood analysis of the transformed COI nucleotide ma¬ 
trix. Numerals are bootstrap percentages for NJ (above branches, 1,000 replicates), MP (below branches, 
1,000 replicates), and ML (parentheses, 100 replicates) analyses. Only bootstrap values greater than 50% 
are shown. 


TABLE 1. Kishino-Hasegawa Test evaluation of the two phylogenetic hypotheses presented in Figure 1 
against the best DNAML tree from the COI nucleotide sequence analysis (Fig. 3). User trees with Z-values 
(= the difference between log likelihood values of best tree and user tree divided by the standard deviation 
of the difference) of absolute magnitude 2.0 or qreater are considered siqnificantly worse (p < 0.05) than the 
best DNAML tree. 


Tree Topology 

Log likelihood value 
(Range of Log 
likelihood values) 

Range of 

differences in Ln L 

Range of Ln L S.D. 

Range of Z-values 

best DNAML 
tree 

-6579.54852 




Figure 1A (34 user 

(- 6661.59438 to 

82.04587 to 

26.3566 to 

-3.1129 to 

trees) 

6634.71091) 

55.16240 

22.8030 

2.4191 * 

Figure 1B (48 user 

(- 6625.17868 to 

-45.63016 to 

27.4627 to 

- 1.6615 to 

trees) 

6601.80601) 

22.25749 

18.5040 

-1.2028 


’Significantly worse than the best tree topology at the p < 0.05 level. 
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FIG. 4. Best tree topology produced by neighbor-joining analysis of the inferred COI amino acid sequences. 
Numerals are bootstrap percentages for NJ (above branches, 1,000 replicates) and MP (below branches, 
1,000 replicates) analyses. Only bootstrap values greater than 50% are shown. 


sels (Unionoida) than to the filibranchous 
mytiloids (mean level of bootstrap support = 
98%). This level of bootstrap support sug¬ 
gests a robust resolution of the evolutionary 
relationships of Neotrigonia among the taxa 
utilized herein (Hillis & Bull, 1993). 

This particular phylogenetic placement of 
Neotrigonia was also manifest in the best 
trees (not shown) constructed from analyses 
of the COI nucleotide sequences when (1) the 
native (untransformed) COI sequences were 
utilized (ML, MP, and NJ analyses), (2) all 
codon positions were coded for transversions 
only (ML, MP, and NJ analyses), and (3) only 
first and second codon positions were used 
(ML, MP, and NJ analyses). These results 
taken together are consonant with the hy¬ 
pothesis that there is a strong phylogenetic 
signal in the COI sequences supporting the 


placement of Neotrigonia as sister taxon to 
the Unionoida. 

Another noteworthy aspect of all of the 
above phylogenetic analyses was the ab¬ 
sence of support for bivalve mollusk mono- 
phyly. In the phylogenetic hypothesis repre¬ 
sented in Figure 3, the protobranch bivalves, 
Nucula and Solemya , are portrayed as a 
clade, with the gastropod. Lepetodrilus, as the 
sister taxon to that clade. The phylogenetic 
hypothesis represented by Figure 4 portrays 
Nucula as the sister taxon of the Autobranchia 
while Solemya is the sister taxon to Lepeto¬ 
drilus . While our analyses were somewhat 
limited due to the relatively small number of 
nucleotides and taxa analyzed, the fact that 
none of the best trees or bootstrap trees gen¬ 
erated from these analyses gave support for a 
monophyletic Bivalvia suggests that the impli- 
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cations of these results be seriously consid¬ 
ered. 


DISCUSSION 

Morton (1987) argued, based partially on 
the anatomical discontinuities between trigo- 
nioids and unionoids, that Neotrigonia was 
a transitional taxon, phylogenetically inter¬ 
mediate between the presumed ancestral 
protobranch bivalves and the more derived 
pteriomorph bivalves. This hypothesis is con¬ 
sistent with that of Salvini-Plawen & Steiner 
(1996; Fig. 1A). However, the results of the 
COI sequence analyses (Figs. 3, 4) strongly 
suggest that the extant representative of the 
Trigonioida, that is, the genus Neotrigonia, is 
the sister taxon to unionoid bivalves, as sug¬ 
gested in the hypothesis of Waller (1990; Fig. 
IB). This phylogenetic propinquity is sup¬ 
ported by similarities in shell structure (Taylor 
et al., 1969, 1973; Tevesz & Carter, 1980), gill 
speculation (Taylor et al., 1969, 1973), sperm 
morphology (Popham, 1979; Healy, 1989), 
and gill cilia patterns (Atkins, 1937; Tevesz, 
1975). Furthermore, the indicated monophyly 
of the Unionoida is consistent with the hy¬ 
pothesis of a single invasion of freshwater by 
the ancestral unionoid. This finding corrobo¬ 
rates the hypothesis of a dramatic evolution¬ 
ary transition in larval morphology, that is, be¬ 
tween glochidium and haustorium/lasidium 
morphology, during unionoid phylogenesis. 
Evaluating the directionality of this character 
state transition will require further, broad- 
scale phylogenetic analyses. 

Another interesting result is the placement 
of Solemya (in Figs. 3, 4) and Nucula (in Fig. 
3), both protobranch bivalves, among the non¬ 
bivalve outgroup taxa. This observation is not 
an artifact of the particular rooting scheme em¬ 
ployed in Figures 3 and 4. It is not possible to 
root either of these topologies such that all of 
the bivalve taxa represented therein form a 
monophyletic group. In all of the phylogenetic 
analyses of the COI sequences, Solemya was 
either (1) the sister taxon to Lepetodrilus or (2) 
in a clade with Nucula and Lepetodrilus . 
Therefore, these analyses provide some sup¬ 
port for the hypothesis that the currently-rec¬ 
ognized molluscan taxon Bivalvia is a poly- 
phyletic assemblage. The non-monophyletic 
status of the Bivalvia was supported by a re¬ 
cent phylogenetic analysis of 18S rDNA se¬ 
quences (Adamkewiczetal., 1997: fig. 2). Fur¬ 
thermore, the topology in our Figure 3 


suggests that the protobranch bivalve genera 
Nucula and Solemya are more closely related 
to the snail, Lepetodrilus, than to the other bi¬ 
valve taxa in the analysis. Thus, the shared 
presence of bipectinate gills and hypobran- 
chial glands (J. E. Morton, 1988), esophageal 
and stomach similarities (Salvini-Plawen, 
1988), ultrafiltration site similarities (Andrews, 
1988), oxygen transport molecule similarities 
(Mangum et al., 1987), and flattened pedal 
areas in both protobranch bivalves and primi¬ 
tive gastropods may be due to shared com¬ 
mon ancestry rather than to the retention of 
ancestral character states. 

The phylogenetic relationships of Neotrigo¬ 
nia, Nucula, and Solemya, as deduced from 
the COI analyses presented herein (e.g., 
Figs. 3, 4), suggest that a significant amount 
of convergent anatomical and conchological 
evolution has taken place within the Mollusca. 
There is a great deal of precedent for this 
statement (e.g.. Allen, 1985; Purchon, 1990; 
Davis, 1994; Salvini-Plawen & Steiner, 1996). 
An important evolutionary implication that 
stems from the phylogenetic placement of 
Neotrigonia as sister taxon to a freshwater 
mussel clade is that it corroborates previous 
hypotheses of convergence (e.g., Waller, 
1990) in the evolution of autobranch bivalve 
gill structure, i.e., either the filibranchous or 
the eulamellibranchous gill condition evolved 
at least twice in the evolutionary history of 
these bivalve taxa (Fig. 5A). The latter possi¬ 
bility is favored by Waller (1990). The evolu¬ 
tion of eulamellibranchous gill organization 
(which facilitated larval brooding) may have 
been a necessary antecedent to the success¬ 
ful colonization and subsequent marked evo¬ 
lutionary diversification in freshwater habitats 
by unionoid bivalves. 

The implications of the hypothesized phylo¬ 
genetic relationships of Nucula and Solemya, 
as inferred from analyses of COI sequences, 
are more profound. It is suggested that the bi- 
valved phenotype has evolved at least three 
times during the evolution of the Mollusca: (1) 
in the ancestor of the Juliidae, a relatively de¬ 
rived family within the opisthobranch gas¬ 
tropods (not represented in the analyses 
herein), (2) in the ancestor of the Autobran- 
chia, and (3) in the ancestor of the Proto- 
branchia (assuming a monophyletic Proto- 
branchia, Fig.5B). If the genera Tuarangia and 
Pseudomyona are found to be bivalved mono- 
placophorans (Runnegar, 1983) rather than 
autobranch bivalves (Mackinnon, 1982; Berg- 
Madsen, 1987), and if the Solemyoida and Nu- 
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FIG. 5. Morphological character optimization using MacClade on the best topology produced from maximum 
likelihood analysis of the transformed COI nucleotide matrix using Katharina as the outgroup. A. The most 
parsimonious estimate of autobranch bivalve gill character state transitions for the taxa included in these 
analyses is that either the filibranchous or eulamellibranchous character state evolved at least twice. B. The 
most parsimonious estimate of character state transitions for the taxa included in these analyses suggests 
that the bivalved phenotype has evolved at least twice. Justification for the use of Katharina as the outgroup 
for the molluscan taxa included in this analysis is provided by numerous studies (e.g., Salvini-Plawen, 1980, 
1985, 1988, 1990: Wingstrand. 1985; Eernisse, et al., 1992; Salvini-Plawen & Steiner, 1996). 
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culoida had independent origins (Purchon, 
1978; Salvini-Plawen & Steiner, 1996), the bi- 
valved condition would have had to evolved at 
least five times within the Mollusca. This strik¬ 
ing assessment is nonetheless not totally un¬ 
expected when evaluated in the context of (1) 
the great evolutionary mutability of body plan 
exemplified by the phylum Mollusca (e.g., J. E. 
Morton & Yonge. 1964; Allen, 1985; Willmer, 
1990) and (2) the multiple origins of bivalved 
external shells in four phyla (Thomas, 1988). 

Under the hypothesis of a polyphyletic Bi- 
valvia, the degree of morphological conver¬ 
gence on the bivalved body plan varies con¬ 
siderably within the extant mollusks. In the 
case of the opisthobranch gastropod genus 
Julia, the shell has converged on the mono- 
myarian (= single adductor muscle) bivalve 
condition while the soft anatomical character¬ 
istics are clearly those of a gastropod (Kay, 
1968). However, the degree of convergence 
in body plans between autobranch and proto¬ 
branch bivalves is much greater and involves 
both shell and anatomical character states. 
The phylogenetic hypothesis displayed in Fig¬ 
ure 3 suggests that in the distinct evolutionary 
histories of the autobranch and protobranch 
mollusks, in addition to the convergent bi¬ 
valved shell, there were convergent (1) losses 
of head and radula, (2)origins of labial palps, 
(3) origins of the dimyarian adductor muscle 
condition, and (4) origins of fibrous ligament. 
Thus, the evolutionary mutability of body 
plans within the Mollusca may be greater than 
that suggested by the traditional classification 
schemes. Phylogenetic evaluations of addi¬ 
tional molecular and morphological data sets 
are needed to test the hypothesized polyphyly 
of the Bivalvia and further decipher patterns of 
molluscan body plan evolution. 
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